The Vlasov model under large magnetic fields in the 

low-Mach number regime 



p. Degond(i'2), S. A. Hirstoaga^^'^'^), M-H. Vignal^^'^) 

May 14th, 2009 

(1) Universite de Toulouse; UPS, INSA, UTl, UTM; Institut de Mathematiques de 

Toulouse; F-31062 Toulouse, France 

(2) CNRS; Institut de Mathematiques de Toulouse UMR 5219; F-31062 Toulouse, France. 

(3) Centre de recherche INRIA Nancy - Grand Est; 615, rue du Jardin Botanique, F-54602 

Villers-Les-Nancy, France. 

Emails: pierre.degond@math.univ-toulouse.fr, sever.hirstoaga@math.univ-toulouse.fr, 

mhvignal@math.univ-toulouse.fr. 

Abstract 

This article is concerned with the kinetic modeling, by means of the Vlasov equation, 
of charged particles under the influence of a strong external electromagnetic field, i.e. 
when e^, the dimensionless cyclotron period, tends to zero. This leads us to split the 
velocity variable in the Vlasov equation into fluid and random components. The latter 
is supposed to have a large magnitude of order 1/e (which corresponds to the low 
Mach number regime). In the limit e — > 0, the resulting model is a hybrid model which 
couples a kinetic description of the microscopic random motion of the particles to a 
fluid description of the macroscopic behavior of the plasma. The microscopic model is 
a first-order partial differential system for the distribution function, which is averaged 
over the ultra-fast Larmor gyration and the fast parallel motion along the magnetic 
field lines. The perpendicular component (with respect to the magnetic field lines) 
of the bulk velocity is governed by the classical relations describing the E x B and 
diamagnetic drifts, while its parallel component satisfies an elliptic equation along the 
magnetic field lines. 
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1 Introduction 



In this article, we are interested in the description of the dynamics of charged particles sub- 
mitted to a large non uniform magnetic field. This problem is of great potential interest for 
describing strongly magnetized plasmas such as those encountered in Tokamak devices like 
ITER. The study of plasma confinement due to a large magnetic field requires to solve the 
Maxwell equations coupled to the description of the plasma turbulent transport. This trans- 
port can be modelled by using either a fluid description [2l[Hlll0l[IIl|16l[nil2ll[2Hll29l[32] 
or a kinetic description |5l [22l [23], [2ll [261 [SHl [M]. Solving three dimensional fluid equations 
is certainly the less expensive way to solve the problem. However, a fluid description usu- 
ally overestimates turbulent fluxes especially in a weakly coUisional regime encountered in 
Tokamaks [9], [18] . Indeed in Tokamaks, the plasma is carried to high temperatures. Thus, 
since the collision frequency decreases with increasing temperature, the plasma enters in a 
nearly collisionless regime. 

By contrast, the kinetic model provides an appropriate description of turbulent transport 
in a fairly general context, but it requires to solve a six dimensional problem (3D in space and 
3D in velocity) which leads to a huge computational cost. To reduce the cost of numerical 
simulations, it is classical to derive asymptotic models with a smaller number of variables 
than the kinetic description (see [3j and references therein). 

Here, we formally derive a new asymptotic model under both assumptions of large mag- 
netic fields and low-Mach numbers. Large magnetic fields usually lead to the so-called drift- 
kinetic limit (see [H [Tj [201 HI] for physics references and [H [6l [131 [El [13 HZ] for mathematical 
results). In this regime, due to the large applied magnetic field, particles are confined along 
the magnetic field lines and their period of rotation around these lines (called the cyclotron 
period) becomes small. However, to our knowledge, the consideration of both large magnetic 
fields and low-Mach numbers is new. As we will see, considering low-Mach numbers brings 
a lot of interesting additional features. 

We consider a simplified plasma model in which we focus on the dynamics of the ions. 
The coupling with the electrons is ignored and the electromagnetic field is assumed to be 
given. In future works, coupling the ion dynamics with those of the electrons and with the 
electromagnetic field is planned. To describe a collisionless ion dynamics (collisions can be 
neglected in Tokamaks in a first instance ; of course, there are situations where collisions 
must be included but we shall discard them in the present work), we use the Vlasov equation. 
In the large magnetic field regime, the Lorentz force term in the Vlasov equation is scaled by 
a large parameter, 1/e^, where stands for the dimensionless ion cyclotron period, i.e. the 
rotation period of the ion about a magnetic field line (or Larmor rotation). The so called 
drift-kinetic or gyro- kinetic regimes are reached when e tends to zero (see [201 [25] ). We shall 
not dwell on the distinction between the drift and gyro kinetic regimes, as we are aiming at 
a different situation. 

Indeed, in addition to the magnetic field being large, we assume that the ion mean 
velocity in the plasma is much smaller than the sound speed. In the Vlasov equation, this 
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assumption implies that the fluid bulk velocity is much smaller than the magnitude of the 
random motion caused by thermal fluctuations [20] ■ Thus, we separate the slow scale of 
the fluid velocity from the fast scale of the random motion and we express the distribution 
function as a function of the random component of the velocity (which is also the velocity 
in the rest frame of the fluid) instead of the velocity in the laboratory frame. In the present 
work, we focus on the case where the Mach number is of order e, i.e. it scales like the square 
root of the dimensionless cyclotron period. This scaling is natural since both quantities 
scale like the square root of the particle mass, and we can view both the large cyclotron 
freqency and small Mach numbers as a consequence of the small particle inertia. It also 
turns out that this scaling hypothesis gives rise to a rich structure in the asymptotic regime. 
The scaled Vlasov equation for the distribution of random velocities must be coupled to the 
fluid momentum equation which provides an equation for the bulk fluid velocity. Therefore, 
the unknowns of the Vlasov model in these new variables are the distribution of random 
velocities and the bulk fluid velocity. The goal of this paper is to investigate the limit e ^ 
of the Vlasov model in this scaling. 

When £ — »• 0, the limit model consists of two sets of equations, one for the distribution 
function of random velocities, and one for the bulk fluid velocity. The distribution func- 
tion of random velocities only depends on space, time and two components of the velocity, 
corresponding to the parallel component along the magnetic fleld line and the magnitude 
of the perpendicular velocity. In other words, the distribution function is independent of 
the (gyro)-phase of the perpendicular velocity in the plane normal to the magnetic fleld 
line. This is a consequence of the ultra-fast cyclotron rotation about the magnetic fleld 
lines. It is convenient to express the distribution of random velocities in terms of the parallel 
velocity and the magnetic moment (or adiabatic invariant), which is proportional to the 
perpendicular energy divided by the magnitude of the magnetic fleld. 

Now, the distribution function in these new variables satisfles a flrst order differential 
system with a constraint. A Lagrange multiplier allows to express this constraint in the 
differential system. The constraint expresses that the distribution function is constant along 
the trajectories of the fast parallel motion along the magnetic fleld lines. This motion is 
characterized by the constancy of the magnetic moment and of some kind of pseudo-energy 
in the parallel direction, which are the two adiabatic invariants of this motion. If a global 
change of variables from the phase space variables to the adiabatic invariants can be found, 
it is possible to eliminate the Lagrange multiplier of the constraint and to express the model 
as a transport equation in the space spanned by these invariants. This transport model 
describes how the electric and pressure forces as well as spatio-temporal variations of the 
magnetic fleld induce a slow evolution of the distribution function function in the space 
of adiabatic invariants. However, it is not always possible to flnd such a global change of 
variables, and, in this situation, the formulation of the problem as a constrained transport 
equation is the only possible expression of the system. Additionally, in most instances, it 
will provide a more flexible formulation for numerical discretization. 

The equation for the bulk fluid velocity is split in two equations, one for the perpendicular 
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component to the magnetic field lines, one for the parallel component. The perpendicular 
component is given through an identity which simply relates this velocity to the E x B- 
drift and diamagnetic drift velocities. The equation for the parallel component is more 
unusual. It is an ellptic equation which expresses how, in the zero-Mach number limit, 
the parallel fluid velocity must adjust in order to guarantee that, at any time, the parallel 
components of the pressure and electric forces balance along the magnetic field lines. This 
elliptic equation is highly anisotropic because posed on each magnetic field line. It is obtained 
by obtained through expressing the constraint of zero parallel force thanks to the moments 
of the distribution function. 

The derivation of the model roughly follows the following steps: we first proceed with 
formal expansions of the two unknowns in powers of the parameter e (Hillbert expansion) 
and we keep the first three orders. The expansion of the momentum conservation equation 
readily leads to an equilibrium constraint expressing that, in the zero-Mach number limit, the 
pressure force must balance the electromagnetic force. This constraint provides an explicit 
relation for the perpendicular component of the bulk fluid velocity. The parallel component 
is only given implicitely through this equilibrium constraint. Finding an explicit equation 
for it requires some moments of the distribution function. The actual computation of this 
equation is slightly involved. 

Now, carrying the Hilbert expansion procedure through for the distribution function 
equation is best done if we change the random velocity variable into a coordinate system 
consisting of the parallel velocity, the energy, and the angle of rotation (or gyrophase) around 
the magnetic field line. Thanks to this coordinate change, we show that the leading order 
term of the distribution function does not depend on the gyrophase. 

Next, we realize that, at each level of the expansion, we are led to inverting the gyrophase 
averaging operator ^T\. We show that the inverse operator can only act on functions 
satisfying a specific solvability condition, namely that their gyrophase average is zero. We 
find the asymptotic model in abstract form by imposing this solvability condition successively 
to the various terms of the expansion, following the classical Hilbert expansion procedure 
of kinetic theory. Providing explicit expressions of the abstract operators appearing as 
outcomes of the Hilbert expansion procedures requires somehow tedious computations, most 
of which will be skipped and given for the reader's convenience in an appendix. 

The remainder of the paper is organized as follows. In section |2l we present the scal- 
ing which expresses the assumptions of strong magnetic field and low Mach number regime. 
Then, we present and comment the main result of this article, namely the asymptotic model. 
In Section [31 by using Hilbert expansions we derive the equilibrium constraint for the lead- 
ing order fluid velocity and we study the equations concerning the leading order distribution 
function. In Section |H we write the abstract asymptotic model and we provide the main 
computational steps which lead to the explicit partial differential system for the limit distri- 
bution function. The conservative form of the model (main result) in terms of the magnetic 
moment is obtained in Section 14.31 Finally, in section 14.41 we give the explicit formula for 
the perpendicular fluid velocity and we obtain the elliptic equation for the parallel part. 



4 



2 The model, the scahng and the main result 



2.1 The Vlasov equation in a strong magnetic field 

We are interested in the dynamics of a single species of positively charged ions in the plasma. 
At this stage of the study, the coupling with the electrons is discarded and the electromag- 
netic field is supposed given. In future work, the model will be expanded by taking into 
account the coupling with the electrons and with a self-consistant electromagnetic field. 

We are interested in finding the asymptotic limit of the Vlasov equation describing the 
dynamics of the ions when they are submitted to a large external magnetic field and where 
additionally, the thermal fluctuations of the velocity are large compared with the bulk fluid 
velocity. 

Denoting by m the ion mass and by q the positive charge of the ion, we start from the 
Vlasov equation 

dtf + v-VJ + ^E + vxB)- VJ = 0, (2.1) 
m 

where / = f{x,v,t) is the distribution function and x ^ fl C. M^, v G M^, and t G M"*" are 
respectively the position, velocity, and time variables. The position x is supposed to belong 
to an open domain f2 C M^. In addition, the electric field E = E{x, t) and the magnetic field 
B = B{x,t) are assumed to be given. 

We supplement this equation with incoming boundary conditions on the boundary dQ 
of fi: 

f{x,v,t) = fB{x,v,t), X e dQ, v-u{x)<0, (2.2) 

where z/(a;) is the outward unit normal to dfl at x and fs is supposed given. represents 
the distribution function of incoming particles in the domain Q. We also prescribe an initial 
datum 

f{x,v,0) = fj{x,v), xen, (2.3) 

where // is the distribution function of particles initially present inside the domain Q. 

Next, we introduce a set of characteristic scales from which an appropriate scaling of 
equation (12. ip will be derived. Let x be a typical length scale of the problem and let v be 
the ion velocity scale (typically v is the thermal velocity of the ions, (2A;gT/m)'^/^, where 
ki3 is the Boltzmann constant and T is the temperature scale). The time scale is therefore 
t = x/v. We denote by 5 G M"*" the characteristic magnitude of the applied magnetic field 
and by £^ = vB that of the electric field. Thus, we define the new variables and given fields 
by 

x' = x/x, v' = v/v, t' = t/t, E'{x',t') = E{x,t)/E, B\x' ,t') = B{x,t)/B. 

Subsequently, letting / the distribution function scale, we introduce the new unknown 
f'ix',v',t') = f{x,v,t)/J. 
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Inserting all these changes into (12.11) and dropping the primes for clarity, we obtain the 
dimensionless equation 

dtf + v- VJ + ^{E + vxB)- VJ = 0. (2.4) 

m 

When the external magnetic field is assumed to be large, the rotation period of the ions 
about the magnetic field lines becomes small. Denoting by ZI^ = ^ the characteristic ion 
cyclotron frequency, we introduce the dimensionless cyclotron period 




Then, under this scaling, the Vlasov equation (12. 4p for f = fe takes the form: 

dtf + V ■ + \{E + vxB)- V.r = 0, (2.5) 
with intial and boundaray conditions still given by (12. 2p and (12.31) . 



2.2 Splitting the drift and fluid velocities 

In many instances, the ion mean velocity in the bulk plasma is much smaller than the sound 
speed [20] . To take into account this observation, we assume that the fiuid ensemble velocity 
is much smaller that the random component of the particle velocity. Therefore, besides 
the magnetic field being large, another key assumption of the present work is that of a low 
Mach number. Because of the large magnetic field, the random component of the velocity 
undergoes a fast motion around the magnetic field lines. We will see that this drives the 
distribution function towards a state which is (at least on the average) isotropic in the plane 
normal to the magnetic field line. This contributes to a reduction of the dimension of the 
problem. 

To implement this idea, we need to decompose the velocity into the fiuid ensemble velocity 
and its kinetic part. For this purpose, we define the local density rf{x,t) and fiuid velocity 
M^(a:, t) as follows: 

7f{x,t)= I f^{x,v,t)dv, n'^u'^{x,t)= I f^{x,v,t)v dv. (2.6) 

The decomposition of the particle velocity is performed through the change of variables 

V = u^{x,t) + c/e, (2.7) 

where c/e is the random component of the particle velocity. Here the random velocity has a 
large magnitude 1/e which is expressed in the choice of the scaling. The distribution function 
expressed in these new variables is denoted by 

r{x,v,t)=e^g'{x,c,t). (2.8) 
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The scaling factor in front of g'^ preserves the density: 

g'^{x,c,t) dc = I f'^{x,v,t)dv. 



Because of the definition of the velocity according to fl2.6p . g^ satisfies the constraint: 

g'{x,c,t)cdc = 0. (2.9) 



Therefore, g'^ does not carry any information about the evolution of the mean velocity itself. 
So, the transformed equation from (12. 5p through the change of variables (12. 7p will consist 
of two parts, an equation for g"^ on the one hand and the momentum conservation equation 
which allows to determine on the other hand. 

After some easy algebraic manipulations, we find that ( 12. 5p is equivalent to the following 
system: 

1 / V ■ 

dtg' + u' ■ V,g' - c ■ {V,u')Vcg' + - c ■ V^g' + ■ V^g' 

e \ 

+ l(c X B) ■ V,/ = 0, (2.10) 
dt{n'u') + V:,. • {n'u' ® u') + ^ (v^ ■ P" - n%E + x 5))= 0, (2.11) 
where VxU'^ stands for the Jacobian matrix {dx^Uj)ij and P^ for the pressure tensor, 

= / g^c®cdc. 



The model is supplemented with incoming boundary conditions: 

/(x,c,t) = 5(|j(x,c,t), xedVt, (u^ + ^c) ■ z/(x) < 0, (2.12) 

where g%{x,c,t) is related to fs through (12. 7p . and with initial conditions 

g'{x,c,Q)=g]{x,c), x e n, (2.13) 

where again, gj is related to // through (12. 7p . 

The statement that eqs. (12. 5p . (12. 6p on the one hand, and (I2.10p . (12. lip on the other 
hand, are equivalent, requires some comment. First, (12. Sp . (12. 6p clearly imply (I2.10p . (12. lip . 
Indeed, eq. (12. lip is simply obtained by multiplying (12.50 by v and integrating over v. Eq. 
(I2.10p is derived by inserting (12. 7p and (12. 8p into (12. 5p and using (12.110 to eliminate dtu'^. 
Conversely, suppose that (I2.10p . (12.110 are satisfied. Then, performing the change of variable 
(12. 7p the reverse way leads to (12.50 . The only thing which remains to be proved is that is 
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the average velocity according to definition (12.61) . This is a consequence that the constraint 
(12. 9p is satisfied by (I2.10p . ( 12. lip . To prove this, let us denote by P{x,t) = J g'^{x,c,t) cdc. 
Then, multiplying (I2.10p by c and integrating over c, we deduce 

e^i^dtF + iu' ■ V.)r + {F ■ V.)m^ + (V. ■ u')F^ = F x B. 

With the initial condition F{x, 0) = 0, we obtain (12. 9p . This shows the claimed equivalence. 

Our goal is to study the asymptotic model formally obtained by taking the limit e ^ in 
(I2.10p . (12. lip . We notice that taking the first moment in (I2.10p and closing with P*^ = p(?T,'^)Id, 
where is the density of 5^'^, Id is the identity matrix and p{n) is a suitable isentropic pressure 
law, we find the isentropic Euler system with Lorentz force, under the scaling used in [8]. 
This remark guided our choice of the present scaling of the Vlasov system. Indeed, it was 
shown in |8] that this scaling allows to derive a sound drift-fiuid model for magnetic plasma 
confinement. 



2.3 Main result: the asymptotic model 

We first introduce some notations. We assume that B does not vanish and we define the 
director of the magnetic field by b = B/\B\. For a particle with random component of the 
velocity c G M^, we introduce its parallel component cy with respect to the magnetic field 
and its magnetic moment /i, given by 

We define also define the parallel and perpendicular components of the bulk velocity u, by 

u = u\\b + u±, u\\ = u ■ b, u± = b X (u X b) . 

The quantity /i is the magnetic moment of the particle in its rotation motion about the 
magnetic field, i.e. the magnetic fiux through the disk enclosed by the circular trajectory of 
the particle. Indeed, (|cp — c|)/2 is the kinetic energy of the transverse motion to B and 
is proportional to the square of the transverse velocity. But the transverse velocity is equal 
to the Larmor radius times the gyrofrequency Uc and r| is proportional to the surface S 
enclosed by the particle motion while Uc is proportional to \B\. Thus 

IcP - d 

' ' I |2 2 2 CI Dl2 

~ |c_l| r^uj^ S \B\ , 

where c± = c — c\\b and ~ means proportionality. Thus, fi ~ S \ B\ which is the magnitude 
of the flux of B through S, as announced. 

For simplicity, we assume that the boundary dQ of Q is a magnetic surface (i.e. at any 
point on the surface, the magnetic field is tangent to it). As a consequence, a magnetic field 
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line starting inside the interior of does not intersect the boundary We also assume that 
the closure f2 of ^2 is a compact set. These assumptions are true for instance in most parts 
of a Tokamak device, since the magnetic field lines are either closed or dense on magnetic 
surface, depending whether the safety factor is rational or not. We recall that a Tokamak 
geometry is that of a torus and each magnetic surface at equilibrium is also a torus. The 
safety factor of a field line is the number of turns around the small circle of the torus which 
are necessary for one turn around the large one, following the field line. If the safety factor 
is rational, the field line is closed otherwise the field line is dense on a magnetic surface. The 
case of magnetic field lines intersecting the boundary can easily be dealt with but will be 
discarded for simplicity. 

The main result of this paper, which will be proven in section H] is the following: 

Theorem 2.1 In the formal limit e ^ 0, any solution {g^,u^) of \2.11\] converges 

to {g,u) given as follows: first, there exists a function G = G{x, fi,c\\,t) such that 

g{x,c,t) = G(x,/i,C||,t). 

Defining Q = 27i\B\G, the functions Q, ux_ and u\\ are solutions to the following problem: 
there exists a function K = K{x, fj,,c\\,t) such that 



^ + 5t^ + ct/c = o, 

at 



c^g:=v.-{ciigb) 

E xb 1 



_d_ 

dc\\ 



((5.$)^) = 0, 



u± 



\B\ 



n\B 



(i.xv.p.) + ?i^f. 

n\B\ 



-3(6- V,)(V, ■ (p||U||6))+2V,. ■ (m||V,. ■ {p\\b)b^ + {E ■ b)V , ■ {nu\\b) + 
+(V, ■ 6) V, ■ ((-3p|| + PA_)u\\b^ +p±(V,. ■ bf u\\ = i?3 , 



(2.14) 
(2.15) 

(2.16) 
(2.17) 



where /C = 2 vr \B\ k. The right-hand side R3, given by (14.371) . does not depend on u\\. Finally, 
the operator S'^ is defined by 



s^g := V. 

d 



dcii 
d_ 



u 



fiVx y< b - b X ^ 



-(V,m) : (6 ® 6) + /iV, ■ f + $ ■ f cng 



d 



— + M ■ VxO In \B\ -Vx-u + V^u ■.{b®b) + — ■ (5 X $) 
at \B\ 



(2.18) 
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with $, f and ¥ given by 

•^ = ^(^--"^-1^1)' f = &x(6-V.)6, w=^^, (2.19) 

\B\ n 

and with the density n and the pressure tensor P defined by 

n = j g{fi,cii)dfxdcii, ¥ = p^{Id-h®h)+p\\b®h, (2.20) 

p_L = J Q{fj,,c\\) fi\B\dfj,dc\\, P\\ = J ^(/^, C||) Cy d/x (ic|[, (2.21) 

and the function K, satisfies 

ICc\\dijdc\\ = 0. (2.22) 



2.4 Comments on the asymptotic model 
2.4.1 The velocity u 

Eq. (12.161) provides an explicit relation for the transverse part of the velocity. In this ex- 
pression, we recognize the E x B drift in the first term, and the diamagnetic drift in the 
second term. The third term is a drift term relating pressure anisotropy and the curvature 
of magnetic field lines. Indeed, (6 ■ Vx)b is proportional to the curvature of the magnetic 
field line times the first normal to the curve. The vector b x {b ■ Vx)b is proportional to the 
curvature times the binormal to the curve (in the Frenet frame). 

Eq. (12.171) is an elliptic equation for the parallel component u\\ of m which is posed along 
the magnetic field line. Indeed, only operators like ip ^ b- Vxf or ip ^ Vx- ipip) are applied 
to u\\ in (12.171) . In the case of closed magnetic field lines, this problem can be uniquely solved. 
In the cae of dense field lines on a magnetic surface, it is a conjecture that this problem can 
be uniquely solved in the space of almost periodic functions. The study of this problem will 
be developped in future work. 

Both equations follow from the relation 

■ P , . 
+ E + uxB = Q. 2.23 

n 

which expresses that, in the zero Mach limit, the total force must vanish. When resolving 
this equation in terms of m^, we immediatly get (12.161) . However, the projection of this 
equation onto the direction parallel to b gives rise to an implicit constraint on u\\: 

V ■ P 

' " •-E)-b = {). (2.24) 



n 



To resolve this constraint into an explicit equation for u\\, one must take derivatives of this 
equation with respect to t and eliminate the time derivatives of ra, and p± which appear 
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from this operation by using velocity moments of the distribution function equation fl2.14p . 
This computation will be detailed at section (I4.4p . Note that a similar procedure was used 
at the level of the drift-fluid limit in 

2.4.2 Trajectories of the fast parallel motion 

To understand the constrained transport model (I2.14p . (12.15^ . it is useful to introduce the 
negative adjoint operator C to C\ given, for a function F = F{x,^,c\\) by: 

CF:=C||6-V,F+(fi-<l>)|^. (2.25) 

Since V^^ ■ S = 0, a simple computation gives 

C^g = C^{2tt\B\G) = 2tt\B\CG = \B\C{^). (2.26) 

jij I 

Therefore, the constraint fl2.15p can be equivalently written 

CG = 0, (2.27) 

and expresses that G = Q/{2tt\B\) is constant on the characteristics of (12.151) or (I2.26p . These 
characteristics are all the curves defined parametrically by [x = X{t), fi = A^(r), cy = G\\{t)) 
where (X(r), A^(r), C||(r)) satisfies the following ODE: 

f^C,,., (2.28) 

dC 

-j^ = b- (f-MVM), (2.29) 
dM 

^ = 0. (2.30) 
dr 

Here, t is frozen and r is a fictitious time which parametrizes the trajectory. The right-hand 
sides of f lOHD - flOOD are evaluated at (X(r), A<(r), q|(r)) and t. 

The integral curves of the ODE system (I2.28p - (l2.30p define the leading order trajectories 
of the particles in their fast parallel motion along the magnetic field lines. Indeed, their 
parallel velocity is G\\ as seen on ( ]2.28p . Similarly, ( 12.30p states that their magnetic moment 
is invariant in this motion. Finally, (I2.29P shows that the parallel force is a combination of 
the projections onto the magnetic field lines of the electric force E on the one hand, and of 
the mirror force —fiVx\B\ on the other hand (see [211 Chapter 4]). This last term takes into 
account the fact that, if the norm of the magnetic field increases along the particle motion, 
more energy is converted into the rotation motion about B, which decreases the magnitude 
of the parallel velocity. Symmetrically, if this norm decreases along the particle path, energy 
is transferred from the rotation to the translation motion and the parallel velocity increases. 
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Now, the constraint (12.271) reflects the fact that the fast parallel motion along the magnetic 
field line is so fast that it instantaneously relaxes G to a constant along these trajectories. 

More precisely, there are three time scales for a particle moving in a large magnetic field: 
the fastest time scale (which we will refer to as the 'ultra-fast' time scale) corresponds to 
the cyclotron or Larmor rotation period about the magnetic field. Then, the second fastest 
scale (referred to as the fast time scale) is the scale of the parallel motion along the magnetic 
field line, and which is described by system (I2.28p - (l2.30p . Then, there is a slow time scale, 
which is that of the various drifts across the magnetic field lines, due to spatio-temporal 
variations of the magnetic field, or to the electric or pressure forces, etc. In this model, 
we focus on the slow time scale and assume that the ultrafast and fast time scales result 
in the homogenization (through time averaging) of the distribution function G. Therefore, 
as a consequence of the avearging over the ultra-fast time scale, G is independent of the 
gyrophase a. And because of the fast parallel motion, G is constant along the trajectories 



While expressing this averaging effect directly on the particle trajectories is difficult, 
the use of the kinetic model directly provides a way to do it by imposing constraints on 
the distribution function. This easier derivation reflects the fact that, to some extent, the 
distribution function describes the particle dynamics in a weak (or equivalently statistical) 
sense. Averaging the trajectories over some fast component is best done by looking at the 
evolution of an observable of the system which is constant over this fast motion. Saying 
that G is independent of the gyrophase and is constant along the trajectories of the fast 
parallel motion is equivalent to saying that there is an equal probability fo find particles 
with different gyrophases, or at different locations along these trajectories. 

2.4.3 Invariants of the fast parallel motion 

The fast parallel motion can be characterized by some invariants. Considering a particular 
magnetic field line, we introduce its curvilinear abscissa s and we write the function i? ■ $ = 
h- {¥ — ixW x\B\) locally along this trajectory as the derivative with respect to s of an effective 
potental V(s, /i): 




(2.31) 



Then, we can write (^^-^^ as 



dX 
ds 



(2.32) 




(2.33) 



dM 



0. 



(2.34) 



ds 
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We immediately deduce two integrals of this motion: 



ci 

M = Constant, + V{s,M) = Constant. 



(2.35) 



The second relation expresses the conservation of a kind of parallel mechanical energy in the 
course of the fast parallel motion. 

Therefore, the phase-space trajectories of the fast parallel motion are supported by mag- 
netic field lines for their space dependence, have constant magnetic moment and are such 
that their parallel velocity is related to the curvilinear abscissa along the magnetic field line 
and to the magnetic moment according to the second relation of (12.351) . The magnetic mo- 
ment is often referred to as the first adiabatic invariant while the second invariant in (12.351) 
is related to the second adiabatic invariant. 

Indeed, the second adiabatic invariant is introduced when the fast parallel motion along 
the magnetic field line has a bounded spatial range. Suppose that we use the second relation 
of (I2.35P to express C\\ as a function of the curvilinear abscissa s and that this relation 
implies that can only be non-negative on a bounded interval [si, S2]. Then, the particle 
is bound to oscillate between the two mirror points Si and S2. This is because, at Si and 
S2, the total energy of the particle, (which is constant along the trajectory) is entirely in 
the form of perpendicular energy. Then, no energy is available for the parallel motion, and 
the only possibility for the particle is to bounce back to the regions where the perpendicular 
energy is lower. This periodic motion is referred to in the physics literature as the bounce 
motion. By classical mechanical calculations, the period of this motion, or bounce period, 
which is also called the second adiabatic invariant, is related to the expression appearing in 
(12.351) . However, this expression has a general meaning, while the second adiabatic invariant 
is only meaningful in the case where the parallel motion is spatially bounded. 

2.4.4 The Lagrange multiplier /C 

Now, /C appearing in (12.141) is the Lagrange multiplier of the constraint (12.151) . Indeed, let 
us consider a test function F{x,fi,c\\) such that 



i.e. such that F is constant along the trajectories of the fast motion. Multiplying (12.141) by 
F and integrating with respect to (x, /i, cy) leads to 



CF = 0, 



(2.36) 



di 



ig,F) + {s^g,F) = o, 




where we have used the notation 




(2.38) 
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Here, Both G and F are constant along the fast trajectories, and it is equivalent to know G 
punctually or to know the collections of all {Q, F) for all test functions F satisfying (I2.36p . 
Indeed, the duality formula (12.381) defines a function EG: 

27r{BG,F) = {g,F), VF such that CF = 0. (2.39) 

Like G, the function BG is constant on the fast trajectories (note that Q = 27t\B\G is not !). 
Now, by the same considerations, {S'^Q, F) determines a function iS^G. This operator is 
unambiguously defined by the relation: 

27r(5tG,F) = (5T^,F), VF such that CF = 0. (2.40) 

Is it difficult to write an explicit expression of the operators BG and S^G without a 
parametrization of the trajectories of the fast motion. As described in section (12.4.31) . this 
requires a parametrization of the magnetic field lines, and a global definition of the two 
invariants (12.351) . Such a global parametrization is not available in general. Section [2.4.51 
below examines the particular case where such a global parametrization exists. 

However, definitions (I2.39P and (I2.40p allow to eliminate the Lagrange multiplier (/C) and 
to write (12.141) according to 

^^{BG)+S^G = 0, (2.41) 

The introduction of the Lagrange multiplier is just a convenient way to express the 
contraint that G is constant along the trajectories of the fast motion, or equivalently, to 
remove the constraint (12.361) on the test function F for the weak form (12.371) of (I2.4ip . 
Indeed, in weak form, system (12.141) . (I2.15P takes the form of a mixed formulation 

j^{g, F) + {S^g, F) + (d/C, F) = 0, VF (2.42) 
{CG,n) = 0, \/n (2.43) 

where now, the two test functions F and Ti are unconstrained. This is the mixed type 
formulation of the constrained problem (12.370 . (I2.36p . We note that /C is defined up to a 
solution of C^JC = 0. It does not seem obvious to single out a more pertinent choice of /C 
among these possible solutions. However, this multiplicity does not affect the problem since 
whatever choice will lead to the same solution Q. 

2.4.5 Expression in the global coordinate system spanned by the invariants (if 
it exists) 

To provide more explicit expressions of the operators BG and S'^G, we assume that there 
exists a global change of variables 

X ^ {y{x),s{x)), (2.44) 
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such that y G M^. A magnetic field fine is defined by y{x) = Constant, and s{x) is the 
curvihnear abscissa along the magnetic field line. We denote by (/i, W) , the two invariants 
defined by fl2.35p and we also assume that W can be defined globally as a function >V(x, /i, cy). 
When the vector-valued function \E'(x, yU, cy) defined by 

^(x, fi, C||) = {y{x),ij, W{x, fi, C||)), (2.45) 

takes a fixed value {yo, fio,Wo), the point cy) in phase-space spans a given trajectory 
of the fast motion denoted by C*(yQ^^o^vVo)- We denote by rf/(?/o,/io,Wo) the length element on 
this curve. 

Now, we obviously can write G as a function of {y, /i, W): 

G{x, /i, C||, t) = Giyix),fi, Wix, /x, cy)). (2.46) 

As a particular test function, we choose 

F{x,fi,c\\,t) = 5(^(x,/x,cy) - (?/o,yUo,>Vo)). (2.47) 

By the coarea formula [12j, the expression involving the delta distribution means that, for 
any smooth (for instance, continuous) test function 0, we have 

0(5(^(a:,/i,cy)-(i/o,/xo,>Vo))rfxrf/it/cy= / (i/o^o,Wo) ^ ^2.48) 



where A = a/ det((i\E' (i\E'*), (i\E' is the derivative of the above defined function \E', d"^* is the 
matrix transpose of \E' and the expression c?\& ci^^* means a matrix product. We note that 
for any function G like fl2.46p . the following identity holds true: 



G(5(^(x,/i,C||) - (?/o, Ato, VVo)) dxd^dc\\ = A(?/o, /io, VVo) ^(i/o, Ato, VVo) (2.49) 
where 

A(?/o,yUo,>Vo) = / 5(^(a;,/i,cy) - (i/o,;Uo,>Vo))rfa;d/idcy. (2.50) 



is the coarea (colength) of C'(y„_^(,^WQ). The function F given by (12.471) obviously satisfies the 
constraint (12.361) . and therefore, fl2.40p leads to the following expression of S'^G: 

(A5tG)(yo,/io,Wo) = J S^g Si^ix, fi,cii) - {yo, fio,m)) dxdfidcii (2.51) 

Now, let us define by G M^, G R, S\\ G M, the fields such that 

s^Q = V. ■ (s^G) + §^{s,g) + -^{S\\g). (2.52) 

15 



We write 

5t6; = V(.,^,,„)-(5e;), 

with S = {Sx, Sfj,, Then, bu using the Green formula and the chain rule: 



(A5tG)(2/o,/io,>Vo) = J V(3;,;.,c||) ■ {Sg)6{-^{x,n,c\\) - {yo, Ho,^o)) dxdndc\\ 

^5* ■ V(^,^,c||) [5(^'(x,/i,C||) - (?/o,/io,>Vo))] dxdfxdcii 

= -J ^V(y,^,w) W) - (z/cyWo, VVo))] |(y,^,w)=*(x,M,c||)(V(^,^,c||)* ■ 5*) dxd^dc\\ 

= ^iyo,t^o,Wo) ■ (^j ^(V(^,/,,c||)^ ■ 5')5(^(a;,/i,C||) - (yo,/io, Wo))rfx(i/i(ic||^ (2.53) 

For the last identity, we have used that, 

V(y,;,,w) W) - (z/o,/io, Wo))] = -V(yo,^o,Wo) [(^((Z/,/^, >V) - (i/o,/io, Wo))] . 

Finally, we can define the vector fields S.y G M^, G M, Syy G M and S = S^, Sw), 
such that 



Sy(?/o,/io, Wq) = y 27T\B\{V^y ■ S:^)6{'^{x, fi,c\\) - {yo, ^q,Wo)) dx dfi dc\\, (2.54) 

S^(?/o,/io,>Vo) = y 27r|S|S'^5(^(a;,/i,C||) - (?/o,/io, Wo))rfxrf/idc||, (2.55) 
Sw(yo,yUo,Wo) = y"27r|S|(V(^,/,,c||)W-5)5(^'(x,/i,C|j) - (?/o,/io,Wo))rfxd/icic||, (2.56) 
and write: 

S^G = A-V(,,^,w) ■ (SG). (2.57) 
Similar considerations lead to the fact that 

B{x, /i, C||, t) = B{y{x),fi, W{x, /x, cy)), (2.58) 

with 

(AS)(2/o,/io,Wo) = 2n\B\S{<i/{x,fi,cn) ~ {yo,ixo,Wo))dxd^dcii. (2.59) 



Finally, system (12.141) . fl2.14p is equivalent to the following system for G{y, fi,yV,t): 

^^{BG) + A-V,,^,w(SG) = 0. (2.60) 
16 



This system is a classical first order hyperbolic equation for G. 

Note that the existence of a global change of variables like (I2.45P is not always guaran- 
teed. Of course, a local one always exist (using the rectification theorem for vector fields 
for instance), but it may not be globally defined. Even if it is the case, the practical deter- 
mination of this change of variables may be a challenging problem, which makes the use of 
these coordinates very often impossible. If the magnetic field evolves in time, this change 
of coordinates must be updated at any new time step, whch increases the computational 
difficulty of the problem even further. 

Another comment is about the finiteness of the integrals (12.501) . (I2.54p to (12.561) and 
(12.591) . If the magnetic field lines are closed, these integrals are on a compact set, and 
they are finite. If the magnetic field lines are not closed, but are dense on a magnetic 
surface, then, by imposing smoothness conditions on the test functions F, we deduce that 
the integrals are not defined on a particular magnetic field line, but rather on the whole 
magnetic surface. Then, we must change the definitions of the coefficients of the system 
accordingly. But with this change, the integrals become surface integrals on a compact 
manifold again (since we assumed from the very beginning that the closure of the domain 
Vl was compact), and are finite as well. Of course, this requires to shift from one definition 
of the coefficients to another one and is very difficult to realize in practice. The weak mixed 
formulation is more suitable because it allows to automatically shift from line integrals to 
surface integrals when passing from closed magnetic field lines to dens ones. 

2.4.6 Initial and boundary conditions 

Let us first make some comments on the geometric assumptions of the boundary dVt, namely 
that dVL is magnetic surface. It implies that no magnetic field line originating from the 
interior of VL intersects dVL. If this were the case, then constaint (12.271) would impose that 
the value of G inside the domain along this field line would be specified by the boundary 
condition at the point where this magnetic field line interesects the boundary. This would 
result in the fact that the limit model would actually be posed on a subdomain Vt' G Vt 
consisting of the union of all magnetic field lines which do not intersect the boundary. In 
order to avoid any further complicated geometrical discussion, we wilil discard this situation 
and assume that dVL is a magnetic field surface. 

The model (12.141) . (12.151) has to be supplemented with suitable boundary conditions on 
dVt. The case where there exist global invariants (y, /x, W) to the trajectories of the fast 
motion will guide us in determining what are the proper boundary conditions. Indeed, in 
this case, the model reduces to the first order differential system (I2.60p in the invariant space 
(l/,/^,>V). 

We first specify some notations. The domain Vl x ]R+ x R where the variable (x, /i, c\\ ) takes 
its values transforms into a domain O through the transformation to the (y, /i, W) variables 
(12.451) . Its boundary is denoted by dO. The boundary conditions must be prescribed on the 
part of do where the vector field E is incoming, i.e. E ■ A < where A denotes the outward 
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unit normal to dO at the considered point on the boundary. 
Let us first assume that g% given by (I2.12p is such that 

g% ^ Gb{x, ^,c\\,t), X e dil, c-iy{x)<0, (2-61) 

as e 0. We also assume that Gb satisfies CGb = 0, i.e. is constant alont the phase- 
space trajectories of the fast motion. This assumption is consistant since, by the geometric 
assumption on dfl, magnetic field lines starting from a point on dQ will remain on dQ. 
Then, we note the following identity: 

(V(,,^,w)-(SG'),F)o + (G,S-V(,,^,w)i^)o= / GF{J: ■ X) dJ:{y, fi,W), (2.62) 

JdO 

where (F, G)o stands for the integral of FG over O and dT^{y, fj,, W) is the superficial measure 
on do. Similarly, we have 

(V(,,^,eii) ■ (Sg), F) + {g,S- V(,,^,e||)F) = [ GF{S^ ■ v) dS{x) d^x rfc,,, (2.63) 

where dS{x) is the superficial measure on dfl and we recall the definition fl2.38p . By the 
calculation developed in section I2.4.5[ the first and second terms of the left-hand side of 
(12.621) are respectively equal to the first and second terms of the left-hand side of (12.631) . we 
deduce that the right-hand sides of these two formulas are equal. Splitting the boundary 
r = dQ X ]R+ X M into r+, r_, and Fq, according to whether [Sx ■ z/) is > 0, < or = 0, and 
splitting do into 9(9+, dO., OOq in a similar way according to the sign of (S ■ A), we get 

/ G'F(S-A)dS(y,/i,>V)- / G'F|(S-A)|dS(y,/i,W) = 

JdO+ JdO. 

= j gF{Sx-v)dS{x)d^idc\\- j gF\{Sx-v)\dS{x)dndc\\. (2.64) 

We now make a very strong assumption, without which it is difficult to prescribe sound 
boundary conditions for this model. We assume that the sign of {S^ ■ i^) is constant along 
the trajectories of the fast motion. Then, it is consistant to assume that the test function F 
satisfies the constraint f l2.36p . vanishes identically on r+ but does not vanish and is actually 
arbitrary on r_. For such a test function, we deduce that 

[ GF|(S ■ A)|rfS(2/,/x, W) = / gF\{Sx- iy)\dS{x)dixdcii. (2.65) 
Jao. Jr. 

Then, prescribing G on dQ. is equivalent to prescribing the value of the integral at the left- 
hand side of (12.651) and, thanks to this relation, to prescribing the integral at the right-hand 
side of (12351) . 
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Therefore, we prescribe the following boundary condition for the solution Q of fl2.14p . 
(EUD: 

j gF\{S^-v)\dS{x)d^idc\\ = j 2-k\B\GbF\{S^ ■ v)\dS{x) d^idc\\. (2.66) 

for any test function F which satisfies the constraint f l2.36p and vanishes identically on r_|., 
and where the datum is given by fl2.6ip . 

If it is possible to use global invariants as coordinates (y,yU, W), we can use (12.650 with 
a test function F of the form (12.470 with (|/o, /xq, VVo) ^ dO_. This leads to 

(G'|(S-A)|dS)(yo,/^o,Wo) = 

27r|5|GB|(5',-z/)|<5(^(x,/i,C||) -(i/o,/io,Wo))c/^(x)d/idc||,, (2.67) 



which gives an explicit prescription for G on the incoming boundary dO^. 

The prescription (I2.66p . or in the cases where it is possible to use global invariants as 
coordinates, (I2.67P are the prescribed boundary conditions for the model (I2.14p . (12.150 . 

The initial conditions are easier. We suppose that the initial conditions (I2.13P are such 
that there exists G/(x, /i, cy), satisfying the constraint (12.271) such that 

g'i Gi, V(x, c) enx R^, (2.68) 

as £ — i> and we prescribe the boundary condition to be 

^1^=0 = Gj., (2.69) 

2.4.7 Comparison with the literature 

Let us introduce the following classical drift velocities due to the variations of the electro- 
magnetic field [21], Chapter 4]: 

Vcd = 1^^' curvature drift due to the magnetic curvature (6 ■ Vx.)6, 

Vgd = "|^(^ ^ ^x\B\), the gradient drift, 
Ved = rn\{E X b), the electric drift, 

\Id\ 

Vd\\ = n(b ■ (Vx X b))b, a drift parallel to the magnetic field. 



19 



With these definitions, the operator iS^ takes the following form 



S^Q ■■= V, • + Kill + + + Vgd)g] 
d 



+ 



dc\\ 
d_ 
9/i 



V 


X 1 


B\ 




B 





dt\B\ 
\B\ 



{b ■ V,)& ■ Ved - V, ■ Ved 



V 


X 1 


B\ 




B 





We recover the classical expressions of the spatial drift which can be found in the literature 
Chapter 4]. By contrast, the drifts in cy or fi which appear here do not appear in general 
in the literature. This work shows that, at least in the regime described by the proposed 
scaling, such drifts must be included otherwise, the limit model is not consistant with the 
Vlasov equation which we used as a starting point. 



3 The asymptotic limit e ^ 0: preliminaries 
3.1 The Hilbert expansion 



Our goal is to find the asymptotic limit e ^ of (12.101) . (12. lip . We start by assuming that 
and admit Hilbert expansions: 



9o + egi + e g2 + 



and = uq + eui + e^U2 + ■ ■ ■ 



Inserting these expansions in equation (I2.10p . we find 



e term: 
term: 

term: 



{cx B) ■ Vcfi-o 
{cx B) ■ Vcgi 

{CX B) ■ Vcfi'2 



0, 

c ■ V^go + Fo ■ Vcgo 
dgo 
dt 



with 



Fn 



no 



Similarly, in (12. lip , we obtain: 



e term: 
term: 



Fi 



Fo 
Fi 



n 



E + uo X B, 
Ui X B, 



(3.1) 
(3.2) 



+ uo ■ - c ■ ((V^Mo)Vc5'o) + 
+c ■ V,f7i + Fo ■ V^gi + Fi ■ V,(7o- (3.3) 
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We will see that the following orders in the expansion of eqs. f l2.10p and (12.111) are not 
required. 

Next, we introduce some formal notations. We denote by L, T, A, and D the following 
operators: 

Lg = -(c X B) ■ VcQ, (3.4) 

Tg = c- V.,g + Fo ■ Ve^?, (3.5) 
dg 

Ag = — + uo-Vxg - c - {(yxUo)Vcg), (3.6) 

Dg = ¥^- V,g. (3.7) 

Then, eqs. (EUD-dM]) take the form 

Lgo = 0, (3.8) 
Lgi = Tgo, (3.9) 
Lg2 = Ag^ + Tg^ + Dg^. (3.10) 

Since the leading order term in the expansion involves the operator L, which describes 
the effect of a circular motion around the magnetic field lines at infinite angular velocity, we 
now specifically examine the properties of this operator. 



3.2 Study of L 

We investigate the solutions of = and more generally, in a second step, of those of 
Lg = h for a given h. Since L operates on c only, the x and t variables are omitted. In what 
follows, we will not seek precise statements about the functional spaces. We assume that 
the functions are as regular as needed for the following statements to be correct. We denote 
by b = B/\B\ the director of the magnetic field. We assume that it is always defined in the 
region of interest i.e. that B does not vanish, and that it is as smooth as needed. 

We first introduce some notations. Let ei, 62 be two vectors such that the set {ei, 62, b} 
forms a direct orthonormal basis of (i.e. b = ei x 62). We denote by (01,02,03) the 
coordinates of c in this basis, with C3 = cy = c ■ b. We denote by c_l = c — c\\b the normal 
component of c to b. It has coordinates (ci,C2,0) in this basis and its norm is given by 
\c±\ = {cf + clY^"^ = (2e — c|)^/^ where e = |cp/2 is the energy. 

Now, we introduce a coordinate system for c which is derived from the cylindrical coor- 
dinates. Any c such that 7^ is uniquely defined by the triple (e, c\\, a) with (e, C||) E V 
and a e where 

I? = {(e, cj|) I e > 0, -v^ < C|| < v^}, (3.11) 
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and S"^ is the one-dimensional torus R/(27rZ), such that 

1 



or conversely 



e = ^(c? + c^ + c^), (3.12) 

C|| = C3, (3.13) 
cos a = ci/{cl + clfl'^, sin a = C2/(c? + c^)^/^ (3.14) 



-J)i/2cosa, (3.15) 
-^)^/^sina, (3.16) 



Cl 


= (2e 


C2 


= (2e 


C3 





(3.17) 

The angle a is the gyrophase while |c_l| is proportional to the gyroradius of the particle in 
the magnetic field. 

Let g{c) = g{e, cy, a) be the expression of g in this coordinate system where the function 
g{e, C||, a) is defined on P x S^. Saying that a G means that g is 27r-periodic with respect 
to a. 

Lemma 3.1 The null space ker L of L consists of functions which only depend on the parallel 
component c\\ = c ■ b and on the energy e = |cp/2, i.e. 

Lg = 3G(e, cy) with (e, cy) G V, such that g{c) = ^(10172 , c ■ 6). (3.18) 



Proof: By elementary algebra, we have 

V =^c+— 6 + — -^^^ 
de c?C|| da 2e — c| ' 

and therefore, taking into account that |6 x c| = |c_l|, we obtain 

r9 = |B||. (3.19) 
Hence, Lg = implies that dg/da = 0, i.e. 

9 = 9{e,c\\), 

which proves fl3.18p . ■ 

In order to solve eq. Lg = h, we introduce the gyroaveraging operator 11 defined for 
every function g{c) by 



% (e, ^) = ^ y ,^ 9ic) 6 (^-^ -ej 6{c ■ b - ci) do = ^ j ^ g{e, c|, a) da, 



(3.20) 
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for all (e, cy) G V. Ug is nothing but the mean value of g over the phase a. 

Since our asymptotic model will be obtained by gyroaveraging the system f l3.8p -( l3?T0l) . 
the following properties of 11 will be useful. Their proof is easy and is left to the reader: 

Lemma 3.2 For any function g{x,c,t), we have: 

n(|^)=j^(n9) + |^.A(nM). (3.21) 

\OXi J OXi OXi OC\\ 

We deduce the following properties, which will be useful in the sequel: 
Lemma 3.3 We have: 

UL = 0, (3.24) 

/ekerL ^ / = n/, (3.25) 

and, for any function g, 

V. + Fo— j ■ U{cg) + (V.6) : —U{c ® c (?) + (Fq • b)^Ug, (3.26) 

UAg= (I- + UO-VJ Ug+ ((^ + uo-vA b] ■ ^Il{cg) 



dt - y - yydt - y j g^i 

-VxUq : ^Ii{c ®cg)- ((V,uo)6) ■ T^n(c^) + (V, ■ mo)%, (3.27) 

oe oc\\ 

UDg = Fi ■ g^ll{cg) + (Fi ■ b)-^ng, (3.28) 
where : denotes the contracted product of two tensors. 

Proof: ( ]3.24p is immediately deduced from ( 13.23^ and the fact that {Bxc)-c= {Bxc)-b = 0. 
(13.251) is obvious from the definition f l3.20p and Lemma 13. 1[ The other formulas are simple 
applications of Lemma 13.21 ■ 

Concerning the moments of g, we have the 

Lemma 3.4 For any function g{x,c,t), we have: 

/ gdc= Ilg{e,c\i)2Trdedcii, 



n 



I gc\\dc= I I\.g{e,c\\) c\\2tt dedc\\ 
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Let g be a function lying in ker L for all {x,t), i.e. such that g = G{x, |cp/2, c ■ b,t). We 
have: 



n = I G{e.,c\\)2Tx dedc 
Iv 



P 



(3.29) 

gc0cdc = p±{Id-biS)b) +p\\b0b, (3.30) 
1 



V 



p±= I G{e,c\\) { e — ^cft ] 2'Kdedc\ 



I G'(e,C||)cj27rciedc||. (3.31) 
Jv 



In particular, we have 
V, ■ P = V,p^ + 



h ■ V,(p|j - pa_) + (pii - Pa_){V^ -b) b+ (pii - p^){b ■ V,)6. (3.32) 



Lemma 3.5 Let h he given. Then, equation 

Lg = h, 

admits a solution if and only if h has zero phase-average i.e. 

Uh = 0. 

// the solvability condition ^3. 34\) is satisfied, all solutions of ^3. 33\) are written 



1 



^(e,C||,a) = — / h{e,c\\,ip)dip + K{e,c\\), 



(3.33) 
(3.34) 

(3.35) 



where K is arbitrary. K can be uniquely determined if we impose to g to satisfy the cancel- 
lation condition 

Ilg = 0. (3.36) 

The so-defined unique solution of l\3. 33\) is denoted by g = L~^h and the operator L^^ is 
called the pseudo-inverse of L. We have 

(3.37) 



L-i/i(e,C||,a) = 1^ J T{a,Lp) h{e, Oil, if) d(p, 



where the Green kernel T{a,(p) is given by 



27r' 



if < ip < a, 



— - 1, if a <(p<27r. 
27r 



(3.38) 



L^^h does not depend on the particular choice of the basis vectors (61,62). If B is a smooth 
function of {x,t), L^^g is a smooth function of {x,c,t) in the open set where B does not 
vanish. 



24 



Proof: Acting 11 onto eq. fl3.33p and using fl3.24p shows that (13.341) is a necessary condition. 
Conversely, using ( I3.19p . we see that any solution of (13.33^ is of the form (13.35^ . The only 
thing to prove is that this formula provides a periodic function of a (otherwise it is not 
possible to invert the change to cylindrical coordinates). This is true precisely if and only if 
h satisfies the solvability condition (13.341) . Applying the cancellation condition (13.361) allows 
to uniquely determine K: 

1 Z"^'' ~ 

^(^'^ll) = ~^r^ / ^(e,C||,(^)(27r-(^)ci(^. 
27r|-«| Jo 

Inserting this expression into (13.351) . we find the expression of L~^h given by (I3.37p . (I3.38p . 
We note that F can be extended by periodicity into a function of 99 defined on S^. Now, if we 
add any constant to F, we do not change the result of (13.370 due to the solvability constraint 
(I3.34p . Precisely, it is an easy matter to see that, changing the origin ei from which we 
measure the phase a amounts to add a fixed constant to a, or to add a fixed constant to F. 
This shows that, despite the use of an origin for the phase a in (I3.37p . actually does not 
depend on this choice. As a result, L~^h only depends on h and has the same regularity as 
6, i.e. it is smooth in the domain where B is smooth and non-zero. ■ 



Remark 3.6 The choice of the vectors 61,62 in the proofs above can be arbitrary. In 
particular, they are not required to form a smooth vector field. We note that this choice is 
equivalent to choose an origin for the phase a. By contrast, c\\ and e do not depend on 61, 62. 
We have shown that the pseudo- inverse despite the fact that its analytical expression 
(13.371) looks dependent on the choice of a particular origin for the phase, actually does not 
depend on it. This remark is important because generating a smooth orthogonal basis 61, 62 
in the plane normal to h is always locally possible but may be globally difficult. Additionally, 
such a basis is not always associated to an admissible coordinate system. Indeed, for this to 
be possible, commutation relations are required. This obviously is not the case in general. 



4 The asymptotic model: derivation 
4.1 The asymptotic model in abstract form 

In this section we derive an asymptotic model for the limit of g'^ by formally passing to 
the limit e — > in (I2.10p - (l2.1ip . This model will be deduced by solving the sequence of 
equations appearing in the Hilbert expansion (13. 80 - 03. lOp . First, by a simple application of 
Lemma 13. ![ eq. (13. Sp is easily solved by: 

Proposition 4.1 There exists a function G(a;, e, cy, t) such that 

goix,c,t) = Gix,\c\y2,c-b,t). (4.1) 



25 



Now, the goal is to find the equation satisfied by go or G. For this purpose, we turn 
to (]3.9p . By a simple application of Lemma [3.51 we have: 

Proposition 4.2 Eq. (13.91) admits a solution gi if and only if go satisfies the solvability 
condition 

T^Tgo = 0. 

// this condition is satisfied, there exists a function k G ker L (in other words, there exists a 
function K{x, e, cy, t) and k = K{x, |cp/2 , c ■ b,t)) such that 

gi = L'^Tgo + k. 

Note that, following (I2.9P which has been proven equivalent to (]2.1ip . we have: 

J g^ C|| do = 0. 

Since, by construction, IIL~^ = 0, and j f do = J Uf do for all functions /, we can write: 

J L-^Tg° Oil dc = j UL-^Tg^cii dc = 0. 
Therefore, we deduce that 

k cii dc = 0. 



Finally, again, by a simple application of Lemma [375| the equation satisfied by go appears 
as the solvability condition of fl3.10p . For such a function, we compute 

Proposition 4.3 Eq. fl3.10p admits a solution g2 if and only if go satisfies the equation 

U{Ago + Tg^ + Dgo) =0. 
We now collect the model in the following theorem and discuss its properties. 

Theorem 4.4 The formal limit e — > o/ problem 1^2. 10\} . ^2. 1 1\) leads to g'^ — > go and 

uo, where go and uo are solutions to the following abstract model: there exists a 
function k such that 

go = Llfyfo, k = Uk, (i.e. % and k are independent of a ), (4-2) 
U{Ago + TL-'Tgo + Tk + Dgo) = 0, (4.3) 
Ii{Tgo) = 0, (4.4) 

V P /V P\ 

Fo = ^ - = E + uo^B, Fi= ^ =Mix5, (4.5) 



no \ n / 1 

kciidc = 0, (4.6) 
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where we recall that the tilde means that function is expressed in the coordinate system \3.12\) - 
\3.1Ji\) and where the operators li, A, T, L"^ and D are respectively defined by (I3.20p . (]3.6p . 
(1331) . (13:371) and fl3TD . 

Let us first note that eqs. (14.21) means that and k do not depend of the gyrophase a. 
In this model, go is determined by eq. (14.31) while the unknown function k plays the role of 
the Lagrange multiplier associated to the constraint (14. 4p . We will see in section that Uq 
is fully determined by the first equation of (14.51) . Finally, it is not necessary to determine 
gi and Ui. Indeed, gi appears in the definition of D but in the next section, using (14.51) . we 
prove that HDgQ = 0. 

In the next sections, we make this model explicit in terms of a partial differential system. 



4.2 The asymptotic model: explicit form 

4.2.1 The constraint UTgo = 

In this section, we consider the constraint (14.41) : 

Proposition 4.5 Condition ( [^.^[ ) is equivalent to the following equation for G: 



Proof: We apply (I3.26P with g = go of the form (14.11) . For such a function, an easy 
computation shows that 



U{cg)=q\Gb, 
Il{c®cg) = G 



[e - -c\\)ild -b®b) + cp< 



(4.7) 
(4.8) 



Then, we insert (gZD, (M into (iX^ by noticing that (V^6) : (b^b) = {{b ■ V^)b) -6 = 0, 
since \b\ = 1, and that (V^fe) : Id = V^; ■ b. ■ 



4.2.2 The main equation (14.3!) 

We start with HAgo- 

Lemma 4.6 We have: 
^Ago 



9G _ , „, ,_ , 
— + V, ■ (uoG) - (V, ■ Mo) ^ 



-e + l4)G 



(4.9) 
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Proof: We insert (iTD, dM]) into IK27h and get 



9 \r. 

-((V.Mo)6) • 7^(c||G6) + (V. ■ uo)G. 

OC\\ 



(4.10) 



We first note that 



d_ 
di 



+ uo ■ V J 6 • 6 = 



since \b\ = 1, which imphes that the second term of (14.101) vanishes. Then, noting that 
(VxUq) : Id = Vx ■ Uq and (VxUq) '■ (6cg)6) = {{'VxUo)b) ■ b, the other terms combine into (14. 9p . 



We now denote 71 = L^^Tqq. With this definition, we have gi = 'ji + k. We now compute 
HT'ji. By inspection of (I3.26p . we reahze that we need to compute n(c7i) and n((c® c) 71) 
(we note that because of (I3.36p . 1171 = 0). We first compute n(c7i): 



Lemma 4.7 We have: 



1 1 
n(c7i) = TbT^^- 2"""^ 



(9 \ f)C 

V, + Fo^jG + C|| — (6-V,)6 



(4.11) 



Proof: postponed in the Appendix. 



Lemma 4.8 We have: 
d 



(Vxb) : — n((c®c)7i) 
1 d 



\B\dcn 



(e-lcj)cii {bxiib-Vx)b))- 



V,. + Fo|-lG 

oe 



■ (4.12) 



Proof: postponed in the Appendix. 

By collecting formulas (14. lip and (14.120 and inserting them into (I3.26p . we obtain the 
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Lemma 4.9 We have: 



nT7i = ( V. + Fo|- ) ■ <{ ^(e - ^cjj) b x 



'dej 1 Ifi 
1 9 



V,. + Fo— G + C|| — (6-V,)6 



9e 



9C| 



1^1 9C|| 



e - icj)c|| (6 X ((6 • V.)6)) • (^V. + Fo^^ G 



Since dk/da = 0, the computation of IITA; is the same as for IIT^fo: 
Lemma 4.10 We have: 

IlTk = (v. + Fo^") ■ {c\\Kb) + (V. ■ b)^ ({e - ^cj) k] + (Fq ■ b)^. 



Finally we calculate IlDgo- 
Lemma 4.11 We have: 



IlDgo = 0. 



Proof: Thanks to ([321), and (|47D, we get: 

nDg„ = Fi . — (n(c(,„)) + (Fi ■ ft)— (%„) 
. (F..6)(|(c„G) + Ag). 

But, with (14.51) . Fi ■ 6 = 0. This concludes the proof. 



4.2.3 The explicit form of the asymptotic model 

From now on, for simplicity, we write u and F instead of uq and Fq, respectively. We collect 
all the previous lemmas in the following theorem, which gives the explicit form of the drift 
kinetic model: 



Theorem 4.12 The formal limit e of problem l{2.10\) - l[2.11\) leads to — > go and 
u'^ u. There exists G such that 



goix,c,t) = go{x,e,cii,a,t) = e, cn, t). 
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and G and u are solutions to the following partial differential system: 

1 



— + V, ■ {uG) - (V, ■ u) — 



dt 



de 



d 



-e + l4)G 



+ v.. + F 



d_ 

de 



; r(e Cm) 6 X 

LB r 2 



J__d_ 

W\dQi 



^c|)c|| (b X ((6 ■ V.)fe)) ■ (v. + F^^ G 



+ 



V. + F^) ■ (c|,ir6) + (V. ■ 6)^ ((e - ic|) iT^ + (F ■ 6) 



dK 

dc\\ 



V. + F^] ■ (c|,G6) + (V. ■ ((e - ic|) g) + (F ■ = 0, 



F 



V.-P 



E + ux 



n 



K c\\ dc = 0, 



0, (4.13) 

(4.14) 

(4.15) 
(4.16) 



where K = K{x,e,c\\,t) is completely determined by the constraint fl4.14p under some suit- 
able boundary conditions, and where we recall the expressions of n and Vx ■ F are given in 
03.291) and 03.321) respectively. 

In the following section, we write the model in terms of more relevant variables. 

4.3 The asymptotic model in the magnetic moment variable 

In order to highlight the physical relevance of the model, it is useful to introduce the new 
variable (see [20j ) 



The quantity fi is the magnetic moment of the particle in its rotation motion about the 
magnetic field, see section [2731 

We introduce the change of variables 



G{x, e, C|[, t) = G{x, /i, C||, t), K{x, e, c\\,t) = K{x, fi, c\\,t) . 



We note that 



{e,c\i)ev^{fi,c\\)ev 



(4.17) 
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where V is defined by (13. lip . 

Since B is divergence free, we liave: 



V, ■ 6 = 5 ■ V J — 



1 

151 



-6- V,(ln|5|). 



Furtliermore, tfie following formulas hold true: 

V,G = V,.G-/i— V,(ln|i?|), 

dG dG cii dG 



dc\\ dc\\ 



\B\ dfi' 



dG 


1 dG 


de 




dG 


dG 


'dt ~ 





9/i dt 



From now on, we work in the variables (/i, cy) and we drop the overbars for clarity. From 
these formulas, we find 



Lemma 4.13 In the new variables, the constraint Iji4-H ) ^■5 written 

,\ dG 



CG := C|| b ■ V^G + b- (y- ij.V^\B\ 



J dc\\ 



0. 



(4.18) 



A function G{x, fi,c\\,t) which satisfies the constraint ( [^.i<^ is constant along the curves 
{x = X(r),/i = M{t), Oil = q|(r)) satisfying the ODE system [MB) - [KM) . 



Lemma 4.14 In the new variables, the asymptotic model Ili4.13\ ), Ili4-14\ ) is written 

dG 



dt 



where C is given by i4-^^ ^^'^ S by 



+ SG + CK = 
CG = 0. 



(4.19) 
(4.20) 



SG := |M + /iV^x6-6x<l> + (-i--^)f I -V^G 



dG 



with 



+ 



■iy^u) : (6®6)+/iV^-f + $-f c 



d 



dc\\ 



dt 



-\Tr,+u-VA\n\B\-V^-u + V^u : (6 ® 6) + — ■ (5 x $) 



\B 



B 



■V. ■ f + $ ■ f 



1 V ■ P 

(F-/iV,.|5|), i = bx{b-W^)b, ¥= " 



15 



n 
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and n and P being given by 

n= G{fi,c\\)2n\B\dfidc\\, ¥ = p^{Id-h®h) + p\\h®h, 
Jv 

p± = G'(/i, C||) yu|i?| 27r|i?| d/i (ic||, P\\ = / C||) c| 27r|i?| d/i (ic||. 

Jv Jv 

Additionally, the velocity u satisfies the constraint: 

¥ = E + uxB. (4.22) 

and K satisfies the constraint 

J Kc\\d(xdc\i=0, (4.23) 

Proof: By performing the change of variables (14.171) into (I4.13p . we readily find that 

dG 

SO ■=u- VccG - VxU : (6 ® 6) cii — 

OC\\ 

- ((^ + ^ ■ In \B\+V^-u- : {h ® 6)^ /i^ 



+ ) ■ \iih X ( V:, + \G \ = \^iW^xh-hx^ \ ■ W^G 



Now, we expand the second order operators and find: 

■ Lb X (v, + $1- 

and 



Cn \ / _ „ _\ dG ( cl „ ^ \ dG 



32 





B\ 




\B 





where we have used that 
and that 

V, ■ (6 X $) - $ ■ (^6 X = ^V. ■ (S X $). 

Collecting these various formulas together leads to fl4.2ip . ■ 
We obtain the conservative form of the model given in Theorem 12.11 defining 

g = 2n\B\G, IC = 2ti\B\K. 

Then, we have the 

Lemma 4.15 The unknowns Q and K, satisfy system 1^2.14^ , Ili2.15\) where 5^ is given 
by Ii2.18\) . $, f and ¥ are defined by ^2.19\) . n and P are given by ^2.20\) . ^2.21\) . Finally the 
velocity satisfies the constraint 

¥ = E + UX B. (4.24) 

Additionally, /C satisfies 



The relations between and iS^ on the hand and C and S on the other hand, are clarified 
in section [2.4.2[ 

4.4 Explicit equations for the limit fluid velocity 

To find explicit equations for the velocity, we will need the moments of Q. We introduce the 
general moments M^,, and Km,q of Q and /C, which are defined by (see [211 Chapter 6]) 

M^^g = y ^ c|p /i" rf/i c/c|| , K^^g = J ICcf fJ.UfX dC\\ . 

We first note that 

n = Mo,o , P\\ = M2,o , p±=\B\Mo,i. (4.25) 

Then, we obtain 
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Lemma 4.16 The moment system satisfied by Mm,q and Km,q is as follows: 



d_ 
at 



6 X F 



u — 



151 



M^,g + ( \B\ (V, X — ) - f ) M„,g+i 



\B\ 



m+2,q 



+ 



¥ 



(m — q) (Vxu) : {h®h) — {m — q) —— • f + g I -77- + u ■ Wx In |i? 



B\ 



d 



dt 



+ qVx-u-^Vx-{bx¥) 



M, 



m,q 



+ 

+ 



-mWx ■ f + (wi — g) 



V 


X 1 


B\ 




B 





■i+^yx-{hxVx\B\) 



m,q+l 



LB 



M 



m,+2,(ji 



+V, ■ {Km+i,qh) - m{h ■ F)ir„_i,, + m{h ■ V,|5|)ir„_i,,+i = 0, (4.26) 
with the constraint 

Vx ■ iMm+l,q b) - m{b ■ ¥)Mra-l,q + m{b ■ V. 1 5 1 )M„_i,g+i = 0, (4.27) 

and with the convention that any moment with negative indexes is identically zero. 

The infinite set of equations fl4.26l) - fl4.27l) for tlie positive integers m and q, contains tlie 
same information as the original model fl2.14p - fl2.15l) . 
In particular, the first moments lead to 

Lemma 4.17 (i) The mass conservation equation holds: 

On 

— + ■ (nu) = . 
(a) We suppose that Mi^ = at the boundary of the domain. Then, 



(4.28) 



Q cii dfi dc\\ = , 



(4.29) 



everywhere. 

(Hi) The constraint lji4-27\ ) for m = 1 and q = carries no information: it is redundant with 
the third equation in Ii2.19\) . 

(iv) The pressures satisfy the following equations: 



dt 



u 



bx¥ 
\B\ 



b \ f ' 

\B\ (Vx X — ) - f j M2,l + — M4,o 



+2 



F 



{Vxu) : (6®6)- — -f 

-D 



2|S|(V,- — )M2,i 



+V. ■ (7^3,0 b) - 2(6 ■ ¥)K,^o + 2(6 ■ Vx\B\)K,^i = 0, 



(4.30) 
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and 



dt 



+ V, 



u - ^^"j p±+( \B\ (V, X A) _ f ) \B\Mo,2 + f M2,i 



+ 



F 6 X F' 

-(V,m) : (6 ® 6) + — . f + V, ■ m - V, ■ — 



+ (v. ■ |5|M2,i + \B\V, ■ 6) = 0. 



(4.31) 



Proof: (i) The equation for the density is obtained by letting m = g = in fl4.26p : 

d 



dt 



bx¥\ L^u^ b 



LB 



B 



+ 



But, using (13.321) . we get: 



b ^ b ^ f 



p± f 

\B\ \B[ 
+V. -(7^1,0 6) =0. (4.32) 



Pi.)- 



Therefore, we have: 

b X n¥ 



and this term is canceled by the divergence operator in (I4.32p . With (12.221) . eq. (14.281) 
follows. 

Point (ii) follows from the application of the constraint (I4.27P for m = and q = 0. 
Indeed, we find 

V, ■ (Mi,o b) = 0, 



out of which (I4.29P follows from the assumption on the boundary conditions. 

Point (iii) also follows from the inspection of the constraint (14.271) but with m = 1 and 
q = 0. Indeed, the left-hand side of this equation is 



(4.33) 



V 


X 1 


B\ 




B 





l.h.s. = V^: ■ {p\\ b) — b ■ n¥ + b ■ , ' p± 

But, thanks to IKW) . IK^Th . we get 

b-n¥ = b-(y,-F) = b- V.pii + (pii - p±){V. ■ b) 
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Therefore, (I4.33P is equal to 



l.h.s. =p^{V,-b) + b- ^j^P± = ^ V, -3 = 0, 

by the divergence free constraint on B. Therefore, the constraint f l4.27p for m = 1 and q = 
is redundant with the definition of F. 

(iv) The equations for the pressures follow from the general moment equation fl4.26p and 

Now, we can turn towards the main result of this section, namely the 
Lemma 4.18 Let Q, /C and u satisfy eqs. i2.14\j , H2.15\) . 1^2. 22\] and ( I4.24p , then 



u = u\\b + u±, (4.34) 

with u\\= u - h and Uj_ = b x [u x b) and u\\ and u± are solutions to (12.160 . (12.170 . 

Proof: In (I4.24p . we insert (14.341) . the decomposition of the velocity into its aligned and 
transverse parts. By taking the vector product of (I4.24p with b and using (I3.32p . we find 
that the transverse part of the velocity satisfies (12.160 . The first component is the classical 
E X B drift. The two other components is the expression of the diamagnetic drift when the 
parallel and transverse pressures are different. 

We now turn to the difficult part: the determination of u\\. For it, we only have an 
implicit constraint, given by the projection of (14.240 onto 6, i.e. 

n{E-b)-(b- V,p|| + {p\\ - p^) (V, ■ 6) ) = . (4.35) 

To show how this leads to a well-posed equation for My, we take the time-derivative of 
(14.350 and use the continuity and pressure equations (14.280 . (I4.30p . (14.310 to eliminate the 
time derivatives of n, p\\ and p^. We first get from (I4.35p : 



df- ' dt' ' dt dt 



d db 1 

+^(P|| - P±) (V. ■ b) + (pii - px)(V, . _) j = , 



which can be written 



^|(^ ■ - I • V.Pii - (Pii -P^)(V. ■ 1^) . (4.36) 
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Now, from fl4.28p . we have 



dn 

'at 



-Vo; ■ {nuiib) - Vx ■ {nu±) 



Using the same methodology with eqs. (14.30p . (14.3ip . we find 



where 



and 



dt 

dp± 
dt 



-"^x ■ {p\\u\\b) - 2p\\ b ■ VccU\\ + Ri, 
-3Vx ■ {p\\u\\b) + 2u\\Vx ■ {p\\b) + Ri , 

-V^ ■ + P±b ■ VxUi\ - p±Vx ■ + R2, 

- ■ - ( ■b) + R2, 



Ri 



u± 



bx¥ 



+ |5|(V. X— -f)M2,l + — M4,o 



\B\ 



+2 



F 



\B\ . 



\B\ 



P|l + 2|5|(V.-^)M2,i 



-V, ■ (7^3,0 b) + 2(6 ■ F)i^i,o - 2(6 ■ V,\B\)K^, 



1) 



R2 



bx¥ 



P± + [ \B\ (V. X — ) - f ) \B\Mo,2 + f M2,i 



+ 



F 6 X F' 

{VxU±) : (6 ® 6) - — ■ f - ■ M± + Vx. 



\B\ 



B 



- (^V, ■ ^) \B\M2,i - \B\Vx ■ (Ki,i b). 

Inserting these formulas into fl4.36p . we find fl2.17p with R^ given by 

R3 = -ib ■ Wx)Ri - (E ■ 6)V, ■ (nu^) - (V,. ■ b) (R, - R^) 



9 , ^ , , db ^ 



db. 

PJ(V..^) 



(4.37) 



This is an elliptic equation for u\\ which is invertible provided boundary conditions for u\\ 
are given at the ends of the magnetic field line. 
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5 Appendix 

5.1 Proof of Lemma [47T 



By (I320D, fl37[5l) - fl3T7j) we have in the frame (ei,e2,6): 

n(^70 = ^ 

1 

But, thanks to (13.19^ . we get: 

If _ , . , I f d 



/ (c7i)(e, 


C||, a) (ia, 








- Cy)^/^ cos a 




- c|)^/^ sin a 




ClI 



(5.1) 



27r 



/" cos a 7i (a) da = — ( (sin a) 71(0;) da, 
J§i 271 Jgi da 



Similarly 



= I sin a— ^ (a) da, 

ZTT Jgi da 

= — ■; — r — / sin a Tfifn (a) da . 

— / sin a 7i (a) da = ■— — / cos a T(yfo (a) da. 
2vr Jgi |i^|27rygi 



27r 

We introduce the new averaging operators, for an arbitrary function h{c) and an arbitrary 
positive integer m: 



W^h{e,c\\) = — / sin(ma) /i(e, cii, a) da, 
2vr Jgi 

1 f 

n^/i(e, C||) = 7^ I cos(ma) /i(e, C|j, a) da. 



which amounts to computing the m-th Fourier coefficients of h with respect to a. Then, the 
previous computation shows that 

/ -(2e - cJ)V2niT^?o 
n(c7i) = p I (2e-c5)V2n^T^o 

The third line corresponds to the applications of the cancellation condition fl3.36p . 

We now need to explicitly compute Tqq (so far, only HTgQ was computed). Using the 
specific form (14. ip of g^, we find: 

Tgo = {¥o-b)— + c- + ¥0— G + (c ® c) : — . (5.2) 



oc\\ V ae / acii 
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Since n^^^(l) = 0, we deduce that 

(d \ dG 
V. + Fo^ j G + n^,5(c ® c) : — , 

and we are left with the task of computing g{c) and g{c® c). Using the same decom- 
position as for (15.11) . we easily find: 

n^c) = i(2e - c||) V2 e^, n^(c) = ^(2e - c\\fl' e„ 

n^(c (g) c) = i(2e - c||)^/2c|| (e2 ® 6 + 6 (g 62), 

n^(c (g c) = i(2e - C||)^/^C|| (ei (g 6 + 6 (g ei). 



Collecting these data, we deduce that: 



and that 



^2e-c5)V^ 

^(2e-4)^/^ 



(V, + Fo9e)G 



(V. + FoSe)^ 



(V.6) 



23 



13 



(V,6)32 
(V.6)31 




-((V.+Fo9e)G)2 
((V. + Fo9e)G)i 





/ -(^-^)23 - (Vx&)32 

+C||^ [ (V,6)i3 + (V,.6)3i 



(5.3) 
(5.4) 



(5.5) 



From the fact that |6| = 1, we have {V xb)h = 0, out of which we deduce that (Vxb)i3 = 0, 
for i = 1,2,3. The first vector in (15. 5p can be easily identified with b x (Vx + ¥ode)G 
while the second one, which reduces to (— (V3;&)32, (Vxfe)3i, 0)"^ (the exponent T denotes the 
transpose), is equal to 6 x ((6 ■ Vx)b). Inserting these last remarks into (15.51) leads to (14.111) 
and ends the proof. 



5.2 Proof of Lemma 14.81 

Using the same method as in the previous lemma, we show that in the basis (ei,e2,&), the 
matrix n((c (g c) 71) has the expression: 



n((c®c) 7i) 



X 

-,2M-r2' 



X 



m 



-(2e-cJ)n|T(7o 

(2e - cl)UlTgo (2e - cJ)n|Tf7o x 

-4(2e-cJ)i/2c||n^T(?o 4(26 - qY/^Ul^Tgo 



39 



where the symbol x indicates that the matrix is symmetric. So, again, we are left with the 
computation of 11^ s^9o- view of O and the fact that, obviously, He sO-) = si^) = 0) 
we need to compute 11^ g{c® c). The same method as previously applies and leads to 



2e-cl 

T 

2e-c' 



n|(c ®c) = — ^ (ei (g) e2 + 62 ® ei) 



,2 



n^(c (g) c) = — - — - (ei (g) ei - 62 ® 62) 



We deduce that 



2e - Cm F)n 

nlTg, = n|(c ® c) : — = — ((V.6)i2 + (V.6) 

rjr; 2e - Cm 

n^T^o = n^(c ® c) : — = — ((V.6)n - 



21 y 



22 j 




9c|| 4 dc\\ 

With fl5.3p and (15.41) . we deduce that Ii{{c®c) 71) has the expression in the basis (ei, 62, 6): 

{2e-clYdG / ~(^^^)i2 - (Vx&)2i X X 

n((c®c)7i) = (V,6)ii - (V,6)22 (V,6)i2 + (V,6)2i x 

^^11 \ 

X 
+ ' ^,J: " I 

-((V.+Fo9e)G)2 ((V. + Foae)^)! 

X X \ 

+ ^ " ^ I Ox. (5.6) 

-(V.6)32 (V.6)31 / 

Now, we need to evaluate (V^6) : n((c(S)c) 71). It is an easy matter to see that the contracted 
product of (Vx&) with the first matrix of (15. 6p is identically zero, as well as with the third 
one. The contracted product of (V^^) and the second matrix of (15. 6p involves the expression 

-(V.6)31 ((V. + Fo9e)G)^ + (V.6)32 ((V. + Fo9e)G)^ = 

= 6-{((V,. + Fo5e)G') x(6-V,.)&}, 

= - ((V. + Fo9e)G) • (fe X (6 ■ V.)6) . (5.7) 

Collecting (15.60 and (15.70 leads to the result (I4.12p and ends the proof of the lemma. 
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